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We study the integer and fractional quantum Hall effect on a honeycomb lattice at half-filling 
(graphene) in the presence of disorder and electron-electron interactions. We show that the inter- 
actions between the delocalized chiral edge states (generated by the magnetic field) and Anderson- 
localized surface states (created by the presence of zig-zag edges) lead to edge reconstruction. As a 
consequence, the point contact tunneling on a graphene edge has a non-universal tunneling exponent, 
and the Hall conductivity is not perfectly quantized in units of /h. We argue that the magneto- 
transport properties of graphene depend strongly on the strength of electron-electron interactions, 
the amount of disorder, and the details of the edges. 
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INTRODUCTION 



Recent progress in development of gate and mag- 
netic field controlled, two-dimensional (2D), graphitic 
devices^'^ has not only opened doors for Carbon micro- 
electronics, but also renewed the interest in the study 
of strongly interacting, low dimensional, electronic sys- 
tems. Graphene is a 2D Carbon material with a hon- 
eycomb lattice and one electron per tt orbital (half-filled 
band), whose elementary excitations are Dirac electrons 
that reside at the corners of the Brillouin zone. These ex- 
citations have linear dispersion relation, e±(k) = ±^'F|k|, 
with a characteristic Dirac- Fermi velocity vp. All elec- 
tronic properties of graphene are determined by the 
physics of Dirac fermions which are quite anomalous 
when compared to the ones found in ordinary electrons: 
the absence of dynamical screening^, a non-Fermi liq- 
uid quasi-particle lifetime^, and anomalous scattering by 
impurities^. Moreover, in the presence of strong disor- 
der, graphite samples (which are obtained from stacking 
of graphene layers) become ferromagnetiofiii indicating 
the important interplay between disorder and electron- 
electron interactions in these materials. We have recently 
shown that because of the low dimensionality, disorder, 
particle-hole asymmetry, and strong Coulomb interac- 
tions, graphene presents the phenomenon of self-doping 
in which extended defects, such as dislocations, discli- 
nations, edges, and micro-cracks, shift the chemical po- 
tential away from the Dirac point to produce electron or 
hole pockets^. The presence of localized disorder, such as 
vacancies and adatoms, leads also to non- trivial physical 
effects that must be understood in order to interpret the 
data correctly. 

In this letter we investigate transport properties of 
graphene under high magnetic fields in the Hall geom- 
etry shown in Fig. ^ When a high magnetic field, 5, 
is applied to a 2D material, the electronic bulk develops 
Landau levels which in the case of Dirac fermions have 
energy = ±i'f^b^V^5 where n is a positive integer, 
Ib = ^y^o/B is the cyclotron radius, and = ch/e 




FIG. 1: (color online) Perspective view of the Hall geometry 
used in our calculations with periodic boundary conditions in 
one direction and zig-zag edges in the other. The magnetic 
field B is applied perpendicular to the graphene plane. 



is flux quanta. Thus, the bulk of the system is gapped 
by the cyclotron energy scale, huOc = V^v^H/Ib (which is 
much larger than the Zeeman energy, qiibB^ where g ^2 
and II B the Bohr magneton^). 

In the integer quantum Hall effect (IQHE), the bulk 
states are gapped and localized due to the disorder, and 
the electronic conduction in a Hall bar occurs through 
its edgesitt. In the case of graphene, due to the Dirac 
fermion nature of its carriers, the Hall conductivity is 
given b}i2*ii: 



^IQHE = {2N + 1 



' h 



(1) 



where N is an integer, and e is the electron charge. Be- 
sides supporting bulk states a graphene Hall bar, such as 
the one shown in Fig.^ also supports surface statesiSii^. 
The result Q is only valid if the surface states do not 
contribute to the conduction. 

In this paper we study the integer and fractional quan- 
tum Hall effect (FQHE) in graphene taking into account 
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the edge and surface states of a graphene Hall bar. We 
show that even in the presence of disorder, when the sur- 
face states become localized, they have a direct effect in 
the magneto-transport. We show that the quantization 
of the conductivity, as given in is not exact by the 
presence of surface states and becomes dependent on the 
details of the sample such as the amount of disorder. Our 
prediction for the FQHE can be verified experimentally 
in graphitic devices^'^. 

The paper is organized as follows: in Section ^ we 
present the Hamiltonian for the problem and show how 
surface and edge Hall modes originate in a graphene Hall 
bar; in Section IHII we discuss the theory for the edge 
modes; Section |^ contains the theory for the surface 
states in the presence of disorder and electron-electron 
interactions; in Section we discuss the Coulomb inter- 
action between edge and surface states and the phase dia- 
gram of a graphene Hall edge in the presence of disorder 
and electron-electron interactions; Section IVII contains 
our conclusions. We have also included one appendix 
with the details of the calculations. 



II. THE MODEL 

The kinetic energy of electrons in graphene is described 
by the Hamiltonian (from now on, we use units such that 
h=l = ks): 

^kinetic = ^ cla^3,CT^t' ^ cJ^^C^,^ + /l.C. , (2) 

where q^ct (cj^r) annihilates (creates) electrons at the site 
Hi with spin a (<j =t, i), t and t' are the nearest neighbor 
and next-nearest neighbor hopping energies, respectively. 
At long- wavelengths the electronic dispersion is given by: 



e±(k) ^ 3t'±VF|k| 



9tV|kP 



where 



3ta 



(3) 



(4) 



In the presence of a magnetic field B the electronic 
hopping between sites and Rj is modified via the 
Peierls substitution, that is, we rewrite the hopping as 
tij exp{z(/?^j} with (^ij = 27r A(r) -dr/^o, where A(r) 
is the vector potential (A = V x B). 

The spectra of the problem for a graphene Hall bar in 
the presence of a magnetic field can be calculated exactly 
by solving an eigenvalue equation^^. In Fig. [21 we show 
the spectra and wavefunctions for a graphene Hall bar 
that is periodic along the direction parallel to the edges 
and 300 lattice spacing wide, such as the one shown in 
Fig.n described in the tight-binding limit with t' = 0.2t 
in the presence of a magnetic field such that the magnetic 
fiux per hexagon, <I>, is ^ = 10~^$o- 



One can clearly see two types of states: bulk states 
that become edge states because of the finite size of the 
Hall bar and a surface state localized at the edge of the 
sample. In the presence of f (that breaks particle-hole 
symmetry) the surface mode is dispersive with a charac- 
teristic velocity: 



vs 



t'a. 



(5) 



while edge states have a velocity: 

v^r^ta, (6) 
where a ~ a/3 x 1.42 A is the lattice spacing. 



wave function squared 




FIG. 2: (color online) Electronic spectra (top) and wave- 
functions (bottom) of a graphene Hall bar as function of the 
momentum k parallel to the edge. For each bullet in the top 
panels (going from left to right), we plot in the bottom panel 
the wavefunction squared as function of the distance x to the 
edge. Note that the two top panels represent the same set of 
bands, (a) and (b): Surface states, (c) and (d): bulk Landau 
levels. Energy given in units of t and distance in units of the 
lattice spacing a. 



III. EDGE STATES 

In order to study the edge states it is conceptually 
simple to consider a large graphene droplet instead of 
a Hall bai»^. When a large magnetic field B = Bz, is 
applied perpendicular to the graphene, a persistent cur- 
rent J fiows along the edge of the droplet, being confined 
by an electric field E created by the termination of the 
graphene droplet. The value of the current is given by: 



J da^yZ X E , 



where 



2e2 



h ' 



(7) 
(8) 
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is the Hall conductivity. Here, v = 2t:1\5^ is the fill- 
ing fraction of the droplet, 5 is the 2D electronic density 
(away from half-filling). Notice, therefore, that the elec- 
trons will drift along the edge with velocity v = Ec/B 
and, hence, the 2D electronic density along the edge, 
Pe{x,t) = 5n{x) (where n{x) is the displacement of the 
edge) obeys the equation: 



dtpe{x, t) - Vdxpe{x, t) = , 



(9) 



which describes a chiral motion (classically, Pe{x^t) = 

Pe{x - Vt)). 

The classical problem can be quantized in terms of the 
Fourier components of the density: 



Pk 



= 7zl' 



(10) 



where L is the circumference of the edge, by canonical 
commutation: 

[Pk,P-k'] = -^kSk^k' . (11) 

ZTT 

The Hamiltonian of the edge waves is then simply: 
27TV sr^ 

H = — yzPkP-k^ 12 

k>0 

The edge fermion operator can be constructed from 
the density operators via a bosonic field, ^(x), such that: 



(13) 



The electron operator, '^^{x)^ is given by the Mandelstam 
construction: 



(14) 



that preserves the commutation relations between the 
electron and the density operators: 



(15) 



This result indicates that the operator carries charge e, 
as required. Furthermore, it is easy to show that: 

^e{x)^e{y) = {-iy^^^e{y)^e{x) , (16) 

and, thus, by the Pauli principle we must require that: 

1^=-, (17) 

m 

where m is an odd integer. The constraint (|T7|) . when 
applied to 0, gives the quantization of Hall conductivity. 
Thus, it is clear that this construction can only describe 
Laughlin's main sequence. For more complicated QHE 
sequences one has to use multiple edge states^. In this 
work we focus on the case given in (|T7)) . 
It is also easy to show that: 



and, hence, 

{^l{x)^em oc e^^^(")^^o)> oc — , (19) 



as expected. In terms of path integrals, the action for 
the ID chiral Luttinger liquid action reads: 



edge 



/ 



m 



dt(t)dx(l) - {dx(l){x)y 



(20) 



where ^(x, t) is a bosonic chiral field along the edge at 
position X and time t. 

This construction stresses the robustness of the Hall 
effect: the edge state being chiral in nature, cannot suf- 
fer any backscattering. It is exactly the electron-electron 
backscattering interaction that creates density wave and 
superconducting states, and the impurity backscatter- 
ing interaction that leads to Anderson localization^^. 
Therefore, chiral edge states iarenot infiuenced either 
by electron-electron interactions or disorder. Finally we 
note that the exponent m is determined entirely by the 
bulk of the system and has topological origin^^. 

In the next section we are going to show that for- 
ward scattering interactions between chiral edge states 
and surface states modify the chiral action introduc- 
ing instabilities in the chiral states at finite momenta. If 
this is the case, it is clear that the relation (|17|) (or ([T^ 
will not hold, spoiling the perfect quantization of the Hall 
conductivity. 




FIG. 3: (color online) Perspective view of edge currents, /, 
on a graphene Hall bar. 



(0(x)^(O)) = —vhi{x) + constant. 



(18) 
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IV. SURFACE STATES AND DISORDER 

The graphene surface state is also a ID state that can 
be described by a non-chiral Luttinger hquid actioni^: 



^surf 



1 



27rK 



I {dtO{x)f 



where 6>(x, t) is a bosonic field at the edge, 
u = vs [l-f/{in^vl)f\ 



(21) 



(22) 



is the renormalized Luttinger liquid velocity (/ is the 
electron-electron forward scattering coupling constant), 
and 

K=[{1- f/{2nvs))l{l + //(27r«s))]'/' , (23) 

is the Luttinger parameter that measures the decay of 
the surface correlation functionsife. The surface density, 
ps{x)^ is written as: 



(24) 



Notice that unlike the edge mode, the surface mode is 
sensitive to the electron-electron interactions and also 
disorder. 

While the problem without disorder can be easily stud- 
ied theoretically^^ (see below), in most of this paper we 
focus on the realistic case of a disordered edge, such as 
the one shown in Fig. ^ Missing Carbon atoms at the 
edge of a graphene hall bar effectively "cuts" the elec- 
tronic wavefunction and leads to lateral confinement of 
the electrons. The lateral confinement of the electrons 
leads to the discretization of the surface state energy lev- 
els, as one would have for a particle moving in a ID box. 
Hence, these electronic puddles have a characteristic ex- 
citation energy scale, or gap, of the order of 



A(^) 



vs t'a 



(25) 



where I is the size of the ID domain. This effect is clearly 
seen in scanning tunneling microscopy (STM) studies of 
the surface of graphite^^. 

The presence of disorder at the edge of the graphene 
Hall bar leads to a back-scattering of the surface electron 
states which can be written as: 



5'bs = 



(y{x)e' 



h.c. 



(26) 



where V{x) is the scattering potential (notice that only 
the 2kF component of the disorder potential contributes 
at low energies, and that the forward scattering part of 
the potential is irrelevant). In the presence of disorder 
V{x) is a random variable with probability: 7^[V^(a:)] (x 
exp { — / dx\V{x)\^ /Vq^ so that, after averaging over dis- 
order: [V(x)V*(y)]^ig^^^^r = Vg8{x - y) . Hence, Vq pro- 
vides a measure of the amount of disorder in the system. 



One possible way of dealing with is via a replica- 
trick^^. In this case, one has to add a new term to the 
Luttinger liquid action (|2T|) : 



^BS 



COS 



[V2{0,{x,t)-0j{x,t'))], (27) 



where Oi{x^t) indicates the field 6 in the i^^ replica. No- 
tice that this term is highly non-local because of the 
quenched disorder. In the absence of edge modes, the 
full action plus ((Tf)) can be understood via a renor- 
malization group (KG) calculation assuming the disorder 
to be weak, that is, we define a dimensionless disorder 
strength, D: 



D 



(28) 



and obtain the RG equations in leading order in 5: 



diD 



D 

2"' 

(3 -2K)D, 



(29) 



where ^ = ln(VFo/^) is the RG scale {W is the running 
energy cut-off of the bosons, Wq ~ t' is the bare cut-off). 
It is easy to see that disorder is irrelevant \i K > 3/2, and 
it is relevant if < 3/2, under the RG fiow. Notice that 
from we have < 1 for repulsive interactions (/ > 
0) and therefore the above RG indicates that disorder 
always fiows to strong coupling, D{i ^ oo) oo, and 
strong interactions, K{i oo) 0, as expected. This 
result implies that the RG breaks down at certain scale 
where D{£'^) ^ 1 and the surface states become Anderson 
localized. The localization scale can be estimated from 
(123 by introducing a localization length, ^, so that: = 
ln(Wo/(^^""^) , and from one finds: 



-l/(3-2K*) 



(30) 



where I^o is the bare amount of disorder in the sys- 
tem and iC* ^ i^(£*). Notice that this result indi- 
cates that there is a characteristic energy scale, E\oc, 
associated with the disorder^ which is of the order of 
^ioc(0 ^ ^s/i ^ t' Dy^'^~'^^ \ Direct comparison with 
shows that ^ can be thought as the typical size of the 
electronic puddles at the edge with a gap in the energy 
spectrum given by £^ioc- 

This result indicates that the bosonic correlations at 
larger distances decay exponentially with ^ (the time cor- 
relations are also short ranged with a characteristic time 
scale Tioc ^ 1/^ioc ^/'^s)- In this case, it is reasonable 
to replace the Luttinger liquid action by: 



(31) 



so that: [(^(x,t)^(0,0))]di,^rder ^ ^We--/^ for x » i 
and t > rioc- 
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Hamiltonian: 

Hll = E qvE{b% + bl^b.,) + Y^ V{q){b% + bl^b.,) , 

q>0 q>0 

(35) 

where V{q) is the surface state forward scattering inter- 
action. Let us consider a generic surface-edge interaction 
potential, V{q), and the interaction Hamiltonian: 

H^,LL = Yl VilHi^l + + ^libl, + b,)] . (36) 

q>0 

The Hamiltonian Hll can by diagonalized via a 
Bogoliubov-Valatin transformation^^, leading to new 
bosonic modes Cq and dqi 



Hll = Y.^(^M^^ 

q>0 



(37) 



FIG. 4: (color online) Schematic representation of a graphene 
Hall bar with a disordered zig-zag edge and its electronic pud- 
dles (see ref. 18]). Squares represent missing Carbon atoms, 
the circles' radii represent the amplitude of the localized sur- 
face electronic wavefunction and the formation of electronic 
puddles. The arrows show the direction of the edge currents. 



V. ELECTRON-ELECTRON INTERACTIONS 
AND EDGE RECONSTRUCTION 

Because the edge and the surface states are confined 
to a small region in space they interact with each other 
via a Coulomb interaction. This interaction leads to a 
new term in the problem: 



^coupling 



dxO dx 



(32) 



with E[q) = y q'^v'^ — [V{q)]'^ . Under the transforma- 
tion the interaction Hamiltonian becomes: 



H^,LL = Yl Viq)Hidl + c,) + a,{d, + c\)] , (38) 

g>0 



with V{q) = V{q)[cosh{Xq) - sinh(A<^)], cosh(2Ag) = 
qvF{q)/E{q), and sinh(2Ag) = V{q)/E{q). The Hamilto- 
nian composed by the sum of (p^ . (pTZ)) . and (|l^8|) . is the 
form of the effective Hamiltonian described in ref. 2l|, 
and can be diagonalized by a generalized Bogoliubov- 
Valatin transformationSiiSi. Introducing a spinor field 
^"f" = (aj, c^q^dq) the total Hamiltonian reads: 



if = ^vl/tZ)vl/-^i?(g), 

g>0 g>0 



(39) 



where D in the grand-dynamical moXih^. The Hamil- 
tonian ((5^ has the form (apart from constant terms): 



where. 



2 



(33) 



is the strength of the electron-electron coupling. 



A. Clean surface states 

Let us consider first the case of a clean surface state 
interacting with a chiral edge state. The edge state is de- 
scribed by the chiral Luttinger liquid Hamiltonian given 
by: 



q>0 



(34) 



where aq (aj) annihilates (creates) a chiral boson with 
momentum g, while the surface state is described by the 



H = '^[u;a{q)alaq + ujf3{q)PlPq + ^i{(l)lllq] , (40) 



q>0 



where, after diagonalization, and the new quasiparticles 
operators read: 



vI>t(Tt)-i = (at 



(41) 



where the matrix T as the form considered in ref. j21ll2^ 

The diagonalization of the Hamiltonian (|39|) amounts 
to find the values of the angles ^, (/>, and r] such that the 
matrix T^^ AT has non-zero diagonal elements only. All 
matrix elements of matrix A are given in appendix H 
As in ref. 21], the relation U0a{q) -^uojs^q) — uo^{q) = qvs 
holds. We have solved the eigenvalue problem for two 
different kinds of electron-electron potentials: {i) a con- 
tact potential given by V{x) = Voa5{x); {ii) V{x) = 
Vo exp(— (|x|/a)). The dispersion of the bosonic modes 
is shown in Fig. O Although both potentials are short 
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FIG. 5: (color online) Quasi-particles spectrum for the two 
potentials considered in the text: (i) dashed line; (ii) soUd 
line. The parameters are Vi = 1.1, V2 = 1.2, V3 = 1.3, 
ai 10.1, a2 ^ 10.2, as 10.3, vs 2.2, ve l- 



ranged, case {ii) introduces a momentum scale ks ^ 
1/a where the spectrum deviates significantly from the 
sound-like behavior obtained with potential (i). It is clear 
from these results that although the coupling between the 
edge modes and surface modes alters the dispersion at fi- 
nite wavelengths it does not lead to any instabilities in 
the clean case. As we are going to show in what follows, 
the presence of disorder changes this picture significantly. 



B. Anderson localized surface states 

If we assume the surface mode is localized as described 
by one can trace the surface mode completely from 
the problem. In fact, using ((ST|) and we find that 
the surface mode is pinned by the edge mode: 



0{x,t) ^ -\i^Ku-^dl(l){x,t) 



(42) 



and therefore, ps{x^t) ex d'^pe{x,t). The pinning of the 
surface mode by the edge mode has a rather interest- 
ing physical interpretation: in the presence of electron- 
electron interactions the surface mode is dragged by the 
edge mode in its motion. The dragging described here has 
similarities with the Coulomb drag between coupled clean 
non-chiral Luttinger liquids in quantum wires^^"^^'^^ but 
differs from it in some fundamental ways: (i) the chi- 
ral edge state is a persistent current, is not subject to 
backscattering, and hence can only interact through for- 
ward scattering (small momentum transfer), as described 
in (p^ : (ii) the non-chiral Luttinger liquid is localized by 
impurities and hence electrons are not free to move (un- 
less the applied force by the chiral component is greater 
than a threshold that leads to the depinning of the lo- 
calized state, a situation not considered here) . Hence, 
although there is no macroscopic voltage drop along the 
edge, the force applied by the edge over the surface state 
will lead to microscopic voltage drops (charge accumu- 
lations in the electronic puddles, see Fig. The final 



picture can be summarized in terms of the scattering of 
the edge electrons by the potential created by the surface 
states. 

Since the theory described by ((^ . (p?T|) . and (p^ is 

gaussian, the surface states can be exactly traced out of 
the problem. The effective action for the chiral modes 
then reads: 



S ■ 



m 



I 



uo - v^k (l- g 



where, 
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muvE 



l<^(fc,c^)l'(,43) 



(44) 



is the surface-edge coupling. Notice that the dispersion 
of the chiral modes is given by: 



u(k) = VEk (l- g 



(45) 



and, hence, for /c <C the chiral mode dispersion be- 
comes: uJk ^ v^ik — nk^) where n = and for k^ 
one finds: uok ^ ve[1 — g]k. 

At long wavelengths {k <^ that is, distances 

larger than the localization length, the surface mode does 
not affect the edge mode. The fact that the dispersion at 
short wavelengths can become negative if ^ > 1 indicates 
the existence of an instability (a quantum critical point) 
at finite wavevectors. It is easy to see that the dispersion 
vanishes at k = kc where: 



1, 



(46) 



for ^ > 1. Therefore, the spectrum of the edge mode be- 
comes negative for k > kc indicating that the edge state 
becomes unstable. Notice that while disorder (and hence, 
^) determines the length scale k~^ of the instability, this 
instability only occurs for a value of g above a critical 
value gc = 1 which marks a quantum phase transition 
in the problem. For g < gc the edge mode is stable but 
for g > gc and for any amount of disorder there is an 
instability in the system with characteristic length scale 
given by The phase diagram of the edge mode as a 
function of the Luttinger parameter, K, and surface-edge 
coupling, ^, has, therefore, the structure shown in Fig. El 
In order to numerically estimate the experimental 
value of this instability, let us consider the case of a 
weakly interacting surface state {K ^ 1, u ^ vs) so that: 



4*. 



Hence, there is a critical density Sc such that g = gc 
given by: 



(47) 
3 = 1 

(48) 



so that k^^ ^ 00 at this point even for small amount 
of disorder. Clearly, in the absence of disorder (I^o = 0, 
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see ((Sn)) ) ^ ^ oo we have = at the outset and the 
instabihty cannot occur. Notice that 5c is independent of 
the disorder, depending only on the ratio between kinetic 
to Coulomb energies in the system. Although there is 
uncertainty^^ on the value of t' ^ if we use t' ~ O.lt ~ 0.2 
eV, ^ 16 eV A, one finds 5c ^ 10^^ - 10^^ cm-^ which 
is the order of magnitude of carriers in these materialsiii. 



K 



3/2 



Stable 



■ Unstable ■ 



g 



FIG. 6: (color online) Phase diagram of the problem, in the 
presence of disorder Do / 0, as a function of the Luttinger 
liquid parameter, K, and surface-edge coupling, g. 

As shown in Ref. 28], this instability is an indica- 
tion of a quantum Hall edge reconstruction where new 
low energy modes are generated at the edge. Edge re- 
construction has been proposed to be important for the 
understanding of IQHE-*^^^^^ as well as FQHE^^ in semi- 
conducting devices and for the interpretation of point 
contact tunneling between a Fermi liquid and a quantum 
Hall edge. In fact, the current- volt age characteristics for 
point contact tunneling follows a power law, / oc F", 
where the exponent a, in the absence of reconstruction, 
is supposed to be universal and independent of the de- 
tails of the edge. Nevertheless, recent experiments show a 
different picture^^ . The discrepancy between theory and 
experiment can be assigned to edge reconstruction. We 
expect a similar effect to occur in graphene and graphite. 



VI. CONCLUSIONS 

The dragging of the surface mode by the edge mode has 
also consequences for the magneto-transport. The longi- 
tudinal, pxx, and Hall, pxy, resistivities depend directly 
on whether the electronic states are localized or not. 
When the electronic states are localized and the longitu- 
dinal conductivity, (Jxx, vanishes, one has pxy = I/c^qhe 



and pxx = and therefore a perfect quantization of the 
Hall resistivity in units of h/e'^. 

We have shown that the presence of edge disorder 
affects directly the physics of the edge states allowing 
for the possibility of edge reconstruction. The effect is 
stronger in the FQHE than in the IQHE because of the 
characteristic energy scales in the problem. In the IQHE 
the dominant energy scale is the cyclotron energy, uJc 
When UJc is larger than the broadening of the Landau lev- 
els due to disorder, the IQHE becomes observable. For 
a magnetic field B ^ 6 T we have cOc ^ 1,000 K, which 
is much larger than the cyclotron energy in conventional 
semiconducting Hall bars (which is of order of a few K). 
This rather large cyclotron energy (a result of the Dirac 
dispersion) makes the observation of the IQHE relatively 
simple^-^'^. 

The situation with the FQHE is very different. For 
the FQHE what matters is the bulk energy of interaction 
between the electrons^^ which is of order of /{cqIb) ^ 
[e^ / {'\/2eQVF)]^c ^ ^c/^o where eo is the dielectric con- 
stant of ffraphene (we have used that e^/{\/2vF) ^ 1, 
see ref. [7|). Because of the presence a back gate in the 
experiments, we expect screening to be as strong as in or- 
dinary semiconducting devices where cq ~ 10 and hence 

/{cqIb) <^ The surface mode localization, and the 
formation of electronic puddles, lead to changes in the 
local electrostatic potentials in the Hall bar and affect 
screening, modifying the bulk electron-electron interac- 
tions if the samples are not wide enough (which is the 
case of the current experiments where the samples are of 
order of 10/i m wideiii). Therefore, the bulk states of the 
FQHE will be directly affected by surface state localiza- 
tion leading to a change of the longitudinal conductivity 
in the system. 

In the presence of a surface state, the longitudinal con- 
ductivity can be small but finite (cTxx <C ctqhe) and 
strongly dependent on the amount of disorder at the 
edge. In this case resistivities are given by: 



r XX (c^xx/c^QHe) , 



^QHE 



1 — (cTxx/ctQHe)' 



(49) 



where p^^^ = I/ctxx- Notice that although there is a 
large reduction in the longitudinal resistivity (pxx <^ Pxx 
since dxx <C ctqhe) it is still finite in the "quantum Hall 
regime" . At the same time, the quantization of the Hall 
resistivity is only partial since it will be spoiled by a 
factor (ctxx/cqhe)^ <^ 1. This type of effect has been 
observed in graphite^^ and we expect it to occur in dis- 
ordered graphene samples. 

In summary, we have studied the integer and frac- 
tional quantum Hall effect in graphene taking into ac- 
count edge and surface modes. We show that although 
the surface modes are localized by disorder in the ab- 
sence of a magnetic field, they become delocalized by the 
edge modes that drag the surface modes in their motion 
via electron-electron interactions. Our results indicate 
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that in this case the Hah edge undergoes a reconstruc- 
tion leading to a non-universal point contact exponent 
that depends strongly on the amount of disorder in the 
system. Furthermore, we also show that a perfect Hall 
effect is not possible in disordered graphene samples due 
to the presence of surface modes. Our results show that 
the Hall resistivity is not quantized and can change sig- 
nificantly from sample to sample depending on disorder 
and electron-electron interactions. 
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APPENDIX A: MATRIX ELEMENTS OF A 

The matrix A has the following diagonal and off- 
diagonal elements: 



1) = qvs cos^ (j) cosh^ - 2V{q) cos 0(sin + cos sinh 6) cosh + E{q) (sin^ (j) + cos^ (j) sinh^ O) , (Al) 

A(2, 2) = {qvs cosh^ 6 sin^ (l)^V{q) cosh 6 (sin(20) - 2 sin^ (j) sinh 6) + E{q) (cos^ + sin^ (j) sinh^ O) ^ cosh^ r] 
+ sinh(27^) ((^(g) cos (j) + {qvs + E{q)) cosh sin (j)) sinh - V{q) cosh(26>) sin (j)) 

+ sinh^ T] (^E{q) cosh^ - 2V{q) sinh cosh + qvs sinh^ 0^ , (A2) 
A(3, 3) = qvs cos^ (f) cosh^ 6 - 2V{q) cos 0(sin + cos sinh 6) cosh 6 + E{q) (sin^ + cos^ sinh^ 6) , (A3) 



^4(1, 2) = V{q) cos(20) cosh cosh 6> + smhr]{{{qvs + E{q)) cos cosh ^ — V{q) sin^) sinh^ — V{q) cos 0cosh(26>)) 

+ ^ cosh 7] sm{2(l)) {qvs - SE{q) + {qvs + E{q)) cosh(2(9) - 2V{q) sinh(2(9)) , (A4) 

A{1, 3) = cos(20) cosh 6> sinh ry - ^ sin(20)(g'?;5 - 3E{q) + (g^5 + E{q)) cosh(26>) 

- 2y(g) sinh(2^)) sinh ry + cosh 7y(y(g) (cos 0cosh(2^) + sin sinh ^) — {qvs -\- E{q)) cose/) cosh smb. 0) , 

(A5) 

A(2, 3) = i sinh(2?7) (2{qvs - 5E{q)) cos^ (f> + {qvs + E{q)){cos{2(j)) - 3) cosh(26i) - 4V{q) cosh6isin(20) 
8 V 

- 2 y(g)(cos(2^) - 3) sinh(26>)^ + ^ cosh{2r]) {2V {q) cosh(2l9) sin(/) 

- {qvs + ^(g)) sinh(26>) sin0 - 2V{q) cos sinh 6>) . (A6) 
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